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A computationally efficient method for calculating the maximum conductance 
of disordered networks: Application to 1-dimensional conductors. 
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Random networks of carbon nanotubes and metallic nanowires have shown to be very useful in the production 
of transparent, conducting films. The electronic transport on the film depends considerably on the network 
properties, and on the inter-wire coupling. Here we present a simple, computationally efficient method for 
the calculation of conductance on random nanostructured networks. The method is implemented on metallic 
nanowire networks, which are described within a single-orbital tight binding Hamiltonian, and the conductance 
is calculated with the Kubo formula. We show how the network conductance depends on the average number 
of connections per wire, and on the number of wires connected to the electrodes. We also show the effect 
of the inter- /intra- wire hopping ratio on the conductance through the network. Furthermore, we argue that 
this type of calculation is easily extendable to account for the upper conductivity of realistic films spanned by 
tunneling networks. When compared to experimental measurements, this quantity provides a clear indication 
of how much room is available for improving the film conductivity. 



I. INTRODUCTION 

The current research on flexible, transparent and con- 
ductive thin films is driven by their huge potential as 
transparent electrodes, with many promising applica- 
tions in the construction of low-power foldable elec- 
tronic displays^—. Metallic nanowires (NW) are the 
most prominent materials of choice to promote electronic 
transport in these films. In this case, electrons move be- 
tween several NW forming a disordered tunneling net- 
work which spans the entire film. Fig. [I] displays a 
microscopy image of one such film. The conductivity 
is primarily limited by the junction resistance that ex- 
ists between neighboring NW. In fact, this has been con- 
firmed by several studies reporting enhancements in the 
film conductivity as a result of improvements in the inter- 
wire contacts^. However, in addition to the inter-wire 
resistance, the network morphology also plays a role in 
limiting the film conductivity. This means that no mat- 
ter how much progress is made in lowering the junction 
resistance, there should be a maximum value for the con- 
ductivity, which is regulated by the network topology. 

From a theoretical perspective, the description of elec- 
tronic transport across these films involves the flow of 
charge through a percolating network made of finite-sized 
1-dimensional (ID) conductors connected to source and 
drain electrodes. Computational studies of this nature 
have so far concentrated on analyzing this phenomenon 
on two very different length scales. The first of these 
approaches is based on the ohmic-behavior displayed by 
percolating networks of finite-length conducting rods. In 
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FIG. 1. SEM image of a silver nanowire film, showing the 
disordered complex network. Courtesy of J. N. Coleman. 



this case the rods are represented by length-dependent 
ohmic resistances, connected to each other by inter-wire 
resistors, and the overall network resistivity is obtained 
by solving Poisson's equation^. A somewhat similar 
method has also been used to describe the film conduc- 
tance when the inter- wire resistance is fully dominant. 
In this case, common when the networks are made of 
carbon nanotubes (CNT) for instance, it is possible to 
solve Kirchhoff's equations to obtain the potential drop 
across the network^. Both methods agree qualitatively 
well with experimental measurements. 

It is also possible to approach the problem from an 
atomistic viewpoint, most notably by means of ab-initio 
calculations coupled with non-equilibrium Green func- 
tion methods. The computational cost of these calcula- 
tions is considerably high, and the absence of periodicity 
on the systems severely limits the studies performed with 
this technique. Nonetheless, recent studies have provided 
results for the intra- and inter-wire conductance in the 
case of a single carbon nanotube junction in the presence 
of 2 and N 2 molecules^. The use of such a method for 
obtaining the quantum conductance between individual 
wires is able to shed some light on the quality of the junc- 
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tion resistances but it is unable to reproduce the general 
observed features since it does not include any disorder 
effects. 



The alternative of carrying out a fully atomistic trans- 
port calculation within a heavily disordered environment 
is very desirable but it is currently too computationally 
demanding and calls for a compromise if one wishes to 
combine the two features. One such a way consists of us- 
ing semi-empirical methods for the electronic structure of 
the wires, e.g. tight-binding models, whose parameters 
are obtained from ab-initio calculation s 12 ' 13 . Because 
this method is significantly less demanding, it is possible 
to include disorder effects which take into account ex- 
tensive configurational averaging procedures. This com- 
bination has been used recently to establish the upper- 
value for the conductivity of CNT-based films with the 
intention of guiding the experimental optimization of 
these films«. By comparing the experimentally mea- 
sured conductivities with the corresponding upper value, 
one knows how much room there is for improvement, usu- 
ally attempted by chemical treatments that can reduce 
the inter- wire junction resistance. This valuable piece of 
information has been used, for instance, to establish that 
the best measured values of conductivities for nanotube- 
based films are too close to the upper limit predicted by 
this technique, which suggests that further improvements 
in their transport properties are unlikely. 



The goal of the present work is therefore to present a 
simple way in which the upper value for the conductance 
of heavily disordered networks can be obtained. Since we 
are interested in the best-case scenario in which the elec- 
tronic transmission is at its maximum, we must eliminate 
potential sources of scattering and decoherence such as 
structural imperfections, impurities, and interaction with 
other quasiparticles. In this situation, it is appropriate 
to consider a purely ballistic regime of transport within 
the network, which calls for a quantum description of 
the conductance. Furthermore, our method is supposed 
to be easily transferable to account for a variety of differ- 
ent materials including nanotubes, noble-metal wires and 
graphenei^, to name but a few. With that in mind, we 
describe the transport calculations within the Kubo for- 
malism expressed in terms of single-particle Green func- 
tions (GF). The advantage of using this method is that 
once the "wires" are chosen, the corresponding GF is 
automatically defined, making results easily transferable 
between different cases. Ballistic transport calculations 
in such an idealized environment provide the maximum 
possible conductance for these disordered systems and 
are easily translated into an upper limit for the conduc- 
tivity of films made of similar networks. In particular, 
these limits can be studied in terms of film porosity and 
in terms of the aspect ratio of the components of the film. 



II. MODEL 

We start by defining the system we are considering and 
how we intend to describe its electronic structure. We 
assume a collection of nanosized elements that are dis- 
orderly oriented and distributed forming a complex net- 
work. This structure is modeled by a random grap h 16 ' 17 , 
which is defined by its total number of nodes (N) and 
the average number of connections per node ((a)). The 
graph is generated such that each node has its own num- 
ber of connections, while the distribution of connections 
obeys a Poisson distribution. Once the interconnected 
network has been generated, a number of nodes ((oe)) is 
chosen to connect the network to each of the source and 
drain electrodes. Notice that in order to avoid short- 
circuits, all nodes in the network are explicitly prevented 
from being in contact with both electrodes at any given 
configuration. 

The electronic structure of individual constituents of 
the system can be described by a tight-binding Hamilto- 
nian of the form 

Up = l/MM/MI + E l/M7b<)8, j\, (1) 

where \ft, i) represents an electron on atom i of an el- 
ement labeled by ft, eo is the on-site energy, 70 is the 
hopping integral, and the sum over (i,j) is over nearest- 
neighboring atoms. This general form can describe sev- 
eral types of nanostructures, depending on the choice of 
atomic structure and on the orbital degrees of freedom 
represented by the states Moreover, it is straight- 
forward to generalize this Hamiltonian to the case when 
there is a collection of N such nanostructures through 
a simple sum Ho = ^fa—xUp as long as there is no in- 
teraction between elements, i.e., assuming they are not 
interconnected. It is simple to connect the nanostruc- 
tures if we introduce a potential bridging any two that 
are in contact. This inter-element coupling appears in the 
form of an electronic hopping term 7' in the Hamiltonian 
given by 

Vw = \MWV',3% (2) 

where ft and ft' label the connecting nanostructures, and 
j and j identify the atoms in each nanostructure. Notice 
that ft and ft' must be different while j and j' are allowed 
to have the same value. 

Associated with the Hamiltonian Up, there are the 
retarded (+) and advanced (—) single-particle GF, defined 
as 

g ± = [E ± -n P r\ (3) 

where E ± = E ± %r\, rj being a small positive imaginary 
part added to the electron energy E. It is worth em- 
phasizing that the choice of Hamiltonian automatically 
defines the corresponding GF. So, for the sake of simplic- 
ity, in what follows we present results for wires consisting 
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of monatomic chains with single orbital atoms. Wires of 
this type are known to carry a single quantum of con- 
ductance ( r = 2e 2 /h) but it is very simple to introduce 
orbital degeneracy to describe wires capable of carrying 
more conducting channels. 

In the case of finite-sized ID atomic chains, it is pos- 
sible to derive an analytic expression for its GF— . The 
(i, j)-th matrix element of the GF of a finite monatomic 
chain with L atoms described by the Hamiltonian in 
Eq. (fTJ) is given by 



g ± 



2 L 
L i E± 

m=l J-J 



sin(i^L)sin(i^) 



e + 2 lQ co S (ff T ) 



(4) 



which is very general and can be used to calculate any 
matrix elements of the GF, as needed. Notice that the 
expression above recovers the familiar result of infinite- 
chain density of states (DOS) in the limit L — >• oo. For 
reasons that will become clear later, we also present an 
analytical expression for the surface matrix element of 
the GF for a semi-infinite monatomic chain, which is 
written as 

2|7or 

where for the ± sign, one has to choose the sign which 
provides a positive DOS. 

Note that Eq.© describes the GF of a single isolated 
wire. The GF of an ensemble of N disconnected wires is 
easily represented by a matrix whose diagonal blocks are 
given by the GF defined in Eq.©. In other words, the 
GF of the disconnected network g ± is written as 
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g ± 



(6) 



where each element in the matrix above is actually a 
matrix itself whose dimension is defined by the length of 
the finite chain representing the NW. The 5-block matrix 
appearing in Eq.© is written for illustration purposes 
only. In general, it is a much larger matrix whose number 
of blocks equals the total number of wires N. The effect 
of the interaction between wires is easily computed using 
Dyson's equation, which can be written as 



G = g ± +g ± VG i 



(7) 



where is the network GF, is the GF of the dis- 
connected wires defined in Eq. © and V is the inter- wire 
coupling matrix. In this particular illustrative example 
the matrix V may look something like 
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(8) 



indicating that wire 1 has a single connection to wire 5; 
wire 2 is connected to wires 4 and 5; wire 4 is linked 
only with wire 2 and wire 5 is connected with wires 1 
and 2. Note that in this pictorial representation wire 3 
is not connected to any of its counterparts, a fact that 
fills the third row and third column of the matrix V with 
zero block matrices. In practical terms, the correspond- 
ing row/column in the matrix will not involve any 
1^-dependent terms and is simply given by Q . It is easy 
to conclude that similar behavior arises in terms of the 
intra- wire indices, that is, a connection between atom j 
of wire /? and atom j' of wire /?' will only involve the 
corresponding matrix elements of the g-matrix. There- 
fore, if a given NW is connected to three neighbors at 
intra- wire sites i, j and fc, the matrix elements Q^, gfp 

gf k , gfj^ G i k , gf^ k (and the respective adjoint elements) 
are the only matrix elements appearing in the expression 
for the network GF of Eq.©. This can be an enormous 
simplification when compared with the full diagonaliza- 
tion method. Rather than inverting enormous matrices 
of the order of N x L, this technique allows us to express 
the network GF in terms of inverse operations of much 
smaller matrices, and whose sizes arc defined by the num- 
ber of connections in the system, i.e., matrices of order 
N x (a). This method becomes extremely advantageous 
when the Kubo formula for the conductance is used, as 
we show in the following section. 



III. KUBO FORMALISM 

The conductance of the system can be calculated by 
the Kubo formalism, which is equivalent to the more fa- 
miliar Landauer method* 9 . The Kubo formalism pro- 
vides a simple expression for the conductance by obtain- 
ing the net current that travels across a reference plane 
in the system, often referred to as the cleavage plane, 
placed between two parts of the system. When this cur- 
rent is expressed in terms of a few GF matrix elements 
of the system, the calculation becomes extremely sim- 
ple. In this formalism the system under study can be 
divided into three parts: a central scattering region (C) 
sandwiched by two perfect semi-infinite leads (L,R), as 
illustrated schematically in the top panel of Fig. [5J The 
location of the cleavage plane, represented by the vertical 
line, is arbitrary and can be chosen in the most conve- 
nient location for the case in question. In our setup the 
conductance of the system can be calculated by plac- 
ing the cleavage plane between either one of the leads, 
and the central region. The leads are represented by 
semi-infinite objects, since it is necessary to have parti- 
cle reservoirs where the charge can move to and from. 
When expressed in terms of GF, this method depends 
only on the GF matrix elements located on both sides 
of the cleavage plane, as if the two halves were totally 
disconnected. 

The implementation presented here is a generalization 
of the usual setup described above, in which the L, C 
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FIG. 2. Schematic representation of the system. The scatter- 
ing region C contains the disordered array of finite nanowires. 
Each perfect lead L and R contain several independent semi- 
infinite wires. On the upper panel the vertical line indicates 
the location of the cleavage plane through which the conduc- 
tance is calculated. 



and R regions are actually made of subsets of the entire 
network. The bottom panel of Fig. [5] shows how the 
system can be separated into leads and a scattering re- 
gion. Analogously to the standard case in which charge 
reservoirs are needed, here we can impose that every wire 
touching the electrode regions L and R is replaced by a 
semi-infinite wire, which provides source and drain for 
the moving charges. The semi-infinite wires in the leads 
do not interact with each other, and are allowed to con- 
nect to only one finite conductor in the network. Re- 
stricting semi-infinite conductors to only one connection 
to the network is not a fundamental requirement, but it 
is a safe and direct way of avoiding short-circuits which 
could provide spuriously high values for the conductance. 
It is also important to bear in mind that the finite wires 
connected to the semi-infinite ones can have numerous 
connections themselves. Translating Fig. [5] into words, 
regions L and R contain {(Xe) semi-infinite wires, and 
region C contains a disordered array of N finite wires, 
where each one has a random number of connections, 
with an average number of connections per wire given by 
(a). 

Having placed the cleavage plane between the regions 
C and R, the Kubo formula for the zero-bias conductance 
is written asS^ 

'Ae 2 ' 



T(E) 



5Re{ 



h 
Tr 



GcVcrGbVcr 



VcrGrcVcrGrc 



}■ 



(9) 



where all the GF are calculated at the Fermi energy Ep, 
and G(E) are the causal GF, defined as 



G(E) = ±-[G-(E)-G + (E)] 



(10) 



It is worth stressing once again the advantages of us- 
ing the set of Eqs. dU> to (J5J) to derive the GF of the 



entire network. As such the matrix elements required to 
evaluate the conductance in the equation above are easily 
calculated. In particular, we emphasize the massive re- 
duction introduced in the GF matrix to obtain only the 
relevant elements. Each finite wire can be represented 
by a GF matrix which contains only the elements corre- 
sponding to atoms interacting with other wires (or the 
leads) and the propagators between those atoms. For 
example, a wire with L atoms is described by a GF ma- 
trix with L 2 elements, but if this wire interacts with only 
a other wires in the network, then its GF matrix can 
be reduced to just a 2 elements. The same reduction is 
applied to the GF of the leads. Each semi-infinite con- 
ductor in the leads is represented by a single matrix ele- 
ment, corresponding to the atom which connects to the 
array. Furthermore, with the use of simple analytic ex- 
pressions for the GF matrix elements, the procedure to 
calculate the network propagators is tremendously sim- 
plified. Moreover, because our expressions are written in 
terms of GF, it is straightforward to replace them with 
GF associated with other types of nanostructures. For 
instance, recent calculations have been carried out with 
numerically evaluated GF for carbon nanotubes^, and 
could possibly be also done for graphenei^. Finally, an- 
other advantage worth mentioning is that the simplicity 
with which we treat the system allows us to evaluate re- 
sults for several samples in parallel. This is paramount in 
calculations performed in disordered environments since 
several samples are required to obtain statistically signif- 
icant results. For each data point presented in the next 
section, several independent configurations were gener- 
ated, and all the error bars (calculated from statistical 
errors in the samples) are smaller than symbol sizes. 



IV. RESULTS AND DISCUSSION 

In order to better understand the role of each param- 
eter on the electronic transport across disordered net- 
works, we study the effects of (a), {cue), and inter- wire 
hopping, independently 

We begin by analyzing the effect of the connectiv- 
ity between wires on the conductance of the network. 
If the number of connections per node is too small, i.e. 
1 < {a) < 2, isolated islands tend to form, instead of 
a connected percolating network, and thus the conduc- 
tance of the system is expected to vanish on the thermo- 
dynamic limit (N — > oo). On the other hand, for {a) = 2, 
w 80% of the wires are interconnected, and a percolating 
cluster is always formed for {a) > 2. However, the elec- 
trodes must be connected to this largest cluster in order 
to form a conducting system. In particular, for (a) = 2, 
there is a 20% chance that any electrode wire will be con- 
nected to one of the smaller isolated islands, and will thus 
not contribute to the electronic transport across the net- 
work. For comparison, we state that for {a) = 4 the frac- 
tion of interconnected nodes is 98%, whereas for {a) = 6 
this fraction becomes 99.5%. 
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FIG. 3. (Color online) Conductance at a fixed energy as a 
function of the connectivity of the network. Electronic trans- 
mission increases as more conducting paths become available. 
Legends indicate total number of wires. The length of the 
wires is fixed at L — 1000, and there are {cue) = 6 connec- 
tions to each electrode. Data averaged over 1000 independent 
configurations. 



In Fig. [3] we show the dependence of the conductance 
at a fixed electronic energy as a function of the connec- 
tivity of the network, for networks with three different 
number of components. The data was obtained consider- 
ing wires of fixed length L — 1000 atoms, and the number 
of wires connected to each electrode (semi-infinite wires) 
was fixed at (cxe) =6. The curves emphasize the influ- 
ence of the network connectivity on the conductance, as 
the former is increased from its lower limit (a) = 1, up 
to values well beyond the percolation limit. Each data 
point is an average over 1000 independent configurations, 
and all realizations were considered in the averaging pro- 
cedure. The increase in conductance is a consequence of 
the increase in the number of available conducting paths 
for the electronic transport across the network. However, 
as more of these conducting paths become available, the 
probability of electrons returning to their original elec- 
trode also increases (zero-bias calculation), which brings 
the system to a steady state where the conductance of 
the network saturates. Furthermore, the increment in 
conductance with the number of wires N seen on Fig. [31 
is also caused by the expansion of conducting paths, since 
the total number of unique connections in the network is 
given by {a)N/2. 

Another major factor on the electronic transport of 
random nanowire networks is the number of wires con- 
nected to the electrodes. It is important to remem- 
ber that each wire connected to an electrode provides a 
charge reservoir from where electrons come, and to where 
they finally go. Moreover, each one of these electrode 
wires supports only one quantum of conductance, which 
introduces a bottleneck effect on the system. For exam- 
ple, if {cue) = 1, no matter how many conducting paths 
are present on the network, only one conduction chan- 
nel is available for electrons moving from the network 
towards the drain electrode, and the maximum conduc- 
tance is therefore limited to one quantum of conductance 



FIG. 4. (Color online) Network conductance as a function of 
the number of connections to electrodes, (ob). Inset: slope of 
the linear-fitted lines increases with the network connectivity 
(a), up to a saturation point. Data is for N = 100, and L = 
1000, averaged over at least 200 independent configurations. 



Tq = 2e 2 /h. On the other hand, for any connectivity be- 
yond the percolation limit several conduction paths will 
be available on the network, thus increasing {cxe) will 
certainly improve the conductance of the network. 

On Fig. 2] we present the direct effect of {cxe) on the 
network conductance, for different values of connectiv- 
ity (a). The conductance of the whole network depends 
linearly on the number of connections to semi-infinite 
conductors, with a slope that varies with the connectiv- 
ity of the network. The dashed lines are linear least- 
square fits to the data. Notice that as {cxe) is increased 
the response in conductance is quite considerable, which 
clearly indicates that once a percolating network is in 
place, the system conductance is limited only by {cxe)- 
On the inset, we show the variation of dT/d{cxE) with the 
connectivity (a), obtained from the slope of the linear- 
fitted lines. The slope reaches a maximum value for large 
enough connectivities, when an increase in the number of 
conducting paths available to the electrons does not alter 
the conductance of the network any further. 

The conductance of random nanowire networks is ex- 
pected to be strongly dependent on the probability of 
electrons to jump between neighboring wires, which is 
related to the inter-wire hopping parameter. This inter- 
wire hopping is defined in analogy to the intra-wire 
hopping integral 70, and is denoted by 7'. In the 
tight-binding model, the hopping integral represents the 
strength of the bond between two atoms in a molecule or 
solid, therefore a higher value indicates a stronger cou- 
pling. It is well known that in metallic nanowire/carbon 
nanotube films, the resistance arises mostly from inter- 
wire electron tunneling, and thus, films with better cou- 
pled elements have lower overall resistance. 

The inter-atomic (intra-wire) hopping is the only pa- 
rameter in the tight-binding model considered here, and 
so it is only natural that we choose to define the inter- wire 
hopping in terms of this parameter. That is to say, we an- 
alyze the influence of the ratio 7V70 011 the conductance 
of the network. If the ratio is too small (7V70 << 1) or 



6 




0.5 1 1.5 2 2.5 3 



FIG. 5. (Color online) Network conductance as a function of 
inter- wire/intra- wire coupling ratio. Maximum conductance 
is obtained around 7V70 = 1.15. N = 100, (a) = 8, («e) = 6, 
L — 1000, and 1000 independent configurations. Inset: Con- 
ductance of two semi-infinite chains connected by a hopping 
7'. Further details on the text. 

too large (7V70 >> 1) the mismatch between the cou- 
plings causes an adverse effect on the electronic trans- 
port through the system, lowering the conductance from 
its maximum possible value. In particular, considering 
only two semi-infinite wires connected to each other at 
the surface by a hopping 7' (forming an infinite linear 
chain), we can see on the inset of Fig. [5] that the maxi- 
mum conductance occurs when the inter-/intra-wire cou- 
plings are equal to each other, i.e., 7V70 — 1- This is 
not surprising if we bear in mind that any contrast in the 
Hamiltonian parameters is a source of additional scatter- 
ing. However, on the main panel of Fig. [5] we show the 
results for network conductance as a function of 7V70 
for networks with N = 100, (a) — 8, (cue) = 6, and 
L = 1000. In this case, it can be observed that the maxi- 
mum conductance occurs for ratios slightly larger than 1, 
which results from the fact that finite-sized wires hanging 
loose off the conducting path are capable of causing simi- 
lar scattering events^. Moreover, we have observed that 
the actual maximum value of T and the ratio at which it 
occurs depends slightly on the fixed value chosen for the 
electronic energy. 



V. CONCLUSION 

In summary, we have presented a general computa- 
tional procedure for the calculation of the quantum con- 
ductance in disordered arrays of nanoscale conductors, in 
the ballistic limit. This method was implemented within 
the tight-binding approximation for the electronic struc- 
ture of individual ID conductors. We have shown how 
the conductance of the network increases with the con- 
nectivity, up to a saturation point. The linear depen- 
dence of the conductance with the number of connections 
to electrodes, for a fixed connectivity, was also presented. 
Finally the dependence of the conductance on the inter- 
wire/intra-wire hopping ratio was calculated, and the op- 



timum point identified. 

In order to calculate the maximum conductivity 
of films made of specific nanostructures (NW, CNT, 
graphene, etc), it is necessary to define a geometric model 
which captures the microscopic arrangement in the film. 
The connectivities of the film can be extrated from the 
model, in terms of quantities like the volume fraction (or 
density) of the film, and the aspect ratio of the individ- 
ual components 10 . In following this recipe we have been 
able to calculate the upper bound for the conductivity 
of films spanned by an array of disordered carbon nan- 
otubes. Most remarkably, this value appears to be pretty 
close to the state-of-the-art experimental measurements 
reported so far—. A similar strategy must be followed in 
order to calculate the upper bound for the conductivity 
of films made of other nanostructures. 
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